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Coulomb suppression of the stellar enhancement factor 
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It is commonly assumed that reaction measurements for astrophysics should be preferably per- 
formed in the direction of positive Q value to minimize the impact of the stellar enhancement factor, 
i.e. the difference between the laboratory rate and the actual stellar rate. We show that the stellar 
effects can be minimized in the charged particle channel, even when the reaction Q value is nega- 
tive. As a demonstration, the cross section of the astrophysically relevant 85 Rb(p,n) 85 Sr reaction 
has been measured by activation between 2.16 < E c .m. < 3.96 MeV and the astrophysical reaction 
rate for (p,n) as well as (n,p) is directly inferred from the data. The presented arguments are also 
relevant for other a- and proton-induced reactions in the p and rp processes. Additionally, our 
results confirm a previously derived modification of a global optical proton potential. 

PACS numbers: 26.50.+X Nuclear physics aspects of novae, supernovae, and other explosive environments, 
24.60.Dr Statistical compound-nucleus reactions, 27.50.-t-e 59 < A < 89 



Introduction. Modern nucleosynthesis studies require 
large reaction networks, often including hundreds and 
thousands of nuclei and their respective reactions with 
light particles. Astrophysical reaction rates employed in 
reaction network calculations are determined either di- 
rectly from cross sections or from the rate for the inverse 
reaction by applying detailed balance. The cross sections 
are either known from experiment or predicted by theory. 
Even when a reaction is experimentally accessible, often 
astrophysical rates cannot be directly measured. Excited 
states are thermally populated in an astrophysical plasma 
whereas only reactions on the ground state of the target 
can be investigated in the laboratory. A measure of the 
influence of the excited target states is given by the stel- 
lar enhancement factor / = r s teilar Ag.s. i defined by the 
ratio of the stellar rate to the ground state rate. 

Relevance of the Q value. The enhancement factor / rev 
for the reverse reaction B(b,a)A (defined by having nega- 
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FIG. 1: Reaction Q values for (p,n) and (a,n) reactions with 

/rev <C /forw - 



tive reaction Q value) is usually larger than the enhance- 
ment /forw of the forward reaction A(a,b)B (being the 
one with positive Q value) because more excited states 
are energetically accessible in nucleus B than in nucleus 
A. This is especially pronounced in photodisintegration 
reactions [l|, Therefore, it was assumed so far that 
more astrophysically relevant transitions are neglected 
when experimentally studying a reaction with negative 
Q value. Furthermore, rates employed in reaction net- 
works often use fitted rates for the forward and backward 
reaction which are derived from the same rate for one di- 
rection to ensure compatibility and numerical stability. 
For the latter reason, the fit for the reverse stellar rate 
is derived from the one for the forward rate by applying 
detailed balance. This is the preferred approach because 
the relation between backward and forward rate @ 
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includes the exponential which would enhance fit inaccu- 
racies when starting from the reverse rate with Q < 0. 

Here, we argue that there are cases for which / rcv < 
/forw due to Coulomb suppression of a part of the ener- 
getically allowed transitions. This effect will be most pro- 
nounced in reactions with a charged particle in one and 
a neutral particle in the other channel, e.g. (n,p), but 
it can also appear when the entrance channel and exit 
channel have Coulomb barriers of different height, e.g. 
(p,cv). Transitions from excited states to the same state 
in a compound nucleus are proceeding at smaller rela- 
tive energy and are stronger suppressed by the Coulomb 
barrier. Thus, a prerequisite is that \Q\ is low compared 
to the Coulomb barrier. Using the NON-SMOKER code 
[H, [3] we calculated /forw an d /rev for reactions involv- 
ing light projectiles (nucleons, a) and targets from Ne to 
Pb between the proton and neutron driplines. We find 
more than 1200 reactions exhibiting the suppression ef- 
fect. Figure [1] shows the obtained range of Q values as a 
function of target charge Z. It can be clearly seen that 
larger \Q\ is allowed with increasing Coulomb barrier. 
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TABLE I: 85 Rb(p,n) 85 Sr reaction decay parameters 



Residual 




Half- 


Gamma 


Relative 7-intcnsity 


nucleus 




lifc 


energy [kcV] 


per decay [%] 


85 E g r 


9/2 + 


64.84 ± 0.02 d 


514.01 ± 0.02 


96 ± 4 


85„g r ™ 


1/2" 


67.63 ± 0.04 m 


231.64 ± 0.01 


84.4 ± 0.2 



Although the strengths of the involved transitions also 
depend on spin and parity of the initial and final state, 
Coulomb repulsion dominates the suppression when the 
interaction energy is small. 

Not only theoretically interesting, the Coulomb sup- 
pression effect is also important for experiments because 
it allows to directly determine an astrophysically relevant 
rate by measuring in the direction of suppressed enhance- 
ment factor. The above mentioned complication of fitting 
rates with negative Q values can be circumvented by di- 
rectly applying detailed balance and numerically com- 
puting the rate for the forward reaction. This is possible 
when / rcv 1. Subsequently, fits for both rates can be 
obtained in the standard way. 

Experimental study of 85 Rb(p,n) S5 Sr. As an example 
of the suppression effect and for the derivation of the 
astrophysical rates, we experimentally studied the reac- 
tion 85 Rb(p,n) 85 Sr, having Q = -1.847 MeV. The impor- 
tance of the reaction is manifold. In the last several years 
a number of proton capture cross section measurements 
with relevance for 7 process studies have been carried 
out (see, e.g., [5j and references therein). The 7 process 
was shown to synthesize p nuclides (proton-rich isotopes 
not accessible to the s and r processes) by a series of 
photodisintegrations of stable nuclides in hot layers of 
massive stars [1, 0, [&, .9]. Recently, systematic 7 process 
simulations found not only that photodisintcgration re- 
actions are important but also that (p,n) reactions, and 
in particular 85 Rb(p,n) 85 Sr, strongly influence the final p 
abundances [l(| ■ Additionally, this reaction is well suited 
to test the optical potential used for calculating the inter- 
action between protons and target nuclei. In low-energy 
(p,7) reactions, cross sections can be sensitive to the pro- 
ton as well as the 7 strength in varying, energy-dependent 
proportions. Conversely, in (p,n) reactions the sensitivity 
is highest to a variation in the proton width because the 
neutron width is much larger, thus eliminating the sen- 
sitivity on the width in the exit channel. Compared to 
the well-studied (p,7) reactions, there is only limited ex- 
perimental information available on the low-energy (p,n) 
cross sections in the mass region of the light p nuclei. The 
cross section of the 85 Rb(p,n) 85 Sr reaction was already 
investigated by [ll[ between E c . m . =3.1 and 70.6 MeV. 
However, the accuracy is not sufficient for astrophysical 
applications, mainly because of the large uncertainty of 
the cm. energies. Moreover, there is only one data point 
in the relevant energy region for the 7 process and it 
bears an uncertainty of ± 0.5 MeV in the cm. energy. 

We measured the 85 Rb(p,n) 85 Sr cross sections using 
the activation method. Some aspects of the experiment 
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FIG. 2: Off-line 7 spectrum measured 9 hours after the end 
of the irradiation at 2.4 MeV. The 514 keV peak from the 
85 Rb(p,n) S5 Sr 9 reaction can be well separated from the anni- 
hilation peak. 

are given here, further details can be found in [5]. Ru- 
bidium targets were produced by evaporating natural 
RbCl onto a thin Al foil. The absolute number of tar- 
get atoms and the uniformity were determined by the 
Rutherford Backscattering Method (RBS) using the Nu- 
clear Microbeam facility of ATOMKI p}. A 2 MeV He+ 
beam was focused down to 3x3 pum 2 and was scanned over 
a surface of 75x75 pum 2 at several different positions of 
the target. Two ion-implanted Si detectors of 18 keV 
system resolution and 50 mm 2 area were applied to col- 
lect the backscattering spectra set at 9 = 165° (ft = 
32.5 msrad) and 9 = 135° (ft = 57.3 msrad). The preci- 
sion of the determination of the number of target atoms 
was better than 3%. The thickness was found to be uni- 
form within 1%. For the cross section measurements the 
targets were bombarded with a proton beam of 600 11A 
typical intensity provided by the Van de Graaff and cy- 
clotron accelerators of ATGMKI. The energy range of the 
proton beam between 2 and 4 MeV was covered with 200 
keV steps. To check a possible systematic errors, the E p 
= 2.6 MeV irradiation was carried out with both accel- 
erators and no difference in the cross section was found. 
Each irradiation lasted approximately 8 hours. 

The 85 Rb(p,n) reaction populates the ground ( 85g Sr) 
and isomeric state ( 85m Sr) of the Strontium isotope, ei- 
ther directly or via 7 cascades. 85g Sr decays by electron 
capture to 85 Rb and 85m Sr with internal transition to 
85g Sr and with electron capture and f3 + to 85 Rb. For 
determining the cross section of 85 Rb(p,n) 85g Sr the 514 
keV, for 85 Rb(p,n) 85m Sr the 232 keV 7 line was used. The 
decay parameters of 85 s> m Sr are summarized in Table I. 
With the exception of a 388.5 keV 7-radiation following 
the decay of 8 Sr, no disturbing 7-radiations induced 
on 87 Rb have been observed since all other open reaction 
channels on 87 Rb lead to stable nuclei. For measuring 
the induced 7 activity a calibrated HPGe detector was 
used. The distance between the surface of the detector 
and the target was 10 cm. This way the 7 yield was 
reasonable and the summing effect was well below 1%. 
The yield of the 511 keV annihilation peak (coming from 
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TABLE II: Experimental cross section and S factor values of the 85 Rb(p,n) 85 Sr reaction. 



Ec.m. 

[keV] 


Ground state cross section 
[mbarn] 


Isomeric state cross section 
[mbarn] 


Total cross section 
[mbarn] 


S-factor 
[lifMeVbarn] 


2158 ± 8 


0.050 ± 0.005 


0.008 ± 0.001 


0.058 ± 0.006 


7.13 ± 0.67 


2341 ± 16 


0.185 ± 0.013 


0.039 ± 0.006 


0.224 ± 0.019 


11.22 ± 0.96 


2552 ± 27 


0.422 ± 0.029 


0.147 ± 0.022 


0.569 ± 0.051 


11.35 ± 1.02 


2566 ± 26 


0.420 ± 0.032 


0.162 ± 0.023 


0.582 ± 0.055 


11.01 ± 1.04 


Ztoo ± Zo 


n oort _L ri 00*7 
0.9z0 ± 0.08 I 


non _L no 1 
V.ZoZ ± O.Ool 


L.ZV ± U.lz 


1U.6D ± 1.00 


2963 ± 30 


1.73 ± 0.17 


0.391 ± 0.039 


2.12 ± 0.21 


9.59 ± 0.92 


3156 ± 32 


3.16 ± 0.29 


0.612 ± 0.061 


3.77 ± 0.35 


9.40 ± 0.87 


3342 ± 35 


4.55 ± 0.40 


1.11 ± 0.14 


5.66 ± 0.54 


8.37 ± 0.79 


3551 ± 36 


7.67 ± 0.64 


1.93 ± 0.23 


9.60 ± 0.87 


8.32 ± 0.76 


3754 ± 37 


11.62 ± 0.90 


2.69 ± 0.32 


14.31 ± 1.22 


7.73 ± 0.66 


3952 ± 40 


14.42 ± 1.18 


5.23 ± 0.64 


19.65 ± 1.82 


6.93 ± 0.64 



proton induced reactions on target impurities) was less 
than or comparable to the one of the relevant transition 
at 514 keV, as shown in Fig. O After each irradiation 
the 7 spectra were taken for 12 h. Because of the rel- 
atively long half life of 85 sSr (T 1/2 = 64.84 d) we were 
able to repeat the activity measurement for each tar- 
get after approximately 1 month when the intensity of 
the 511 keV radiation is substantially reduced. The two 
measurements yielded consistent cross sections proving 
the proper separation of the 511 keV and 514 keV peaks. 

The measured total cross sections cover 3 orders of 
magnitude, varying from 0.058 to 19.645 mb. Table II 
lists the measured cross sections a and the S factors, de- 
fined as S(E) = oE~ x exp(— 2ttt]), with the Sommerfeld 
parameter 77 accounting for the Coulomb barrier pene- 
tration The quoted uncertainty in the E c m . values 
corresponds to the energy stability of the proton beam 
and to the uncertainty of the energy loss in the tar- 
get. The uncertainty of the cross section values is the 
quadratic sum of the following partial errors: efficiency 
of the HPGc detector (6%), number of target atoms (3%), 
current measurement (3%), uncertainty of decay param- 
eters (<4%) and counting statistics (0.7-4%). 

Discussion. The measured S factors are compared 
to theoretical predictions obtained with the code NON- 
SMOKER d, 13 in Fig. [3] The standard calculation 
applied a proton optical potential widely used in astro- 
physical applications, based on a microscopic approach 
utilizing a local density approximation [15| ]. Low-energy 
modifications, which are relevant in astrophysics, have 
been provided by [l6[ . As can be seen in Fig. [3l the the- 
oretical energy dependence of the resulting S factor is 
slightly steeper than the data, although there is general 
agreement in magnitude. In the energy range covered by 
the measurement, the proton width is smaller than the 
neutron width (except close to the threshold) and thus 
uncertainties in the description of the proton width (and 
proton transmission coefficient) will fully impact the re- 
sulting S factor. A recent investigation Q suggested that 
the strength of the imaginary part of the microscopic po- 
tential should be increased by 70%. We find that the 
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FIG. 3: Experimental (dots) and theoretical (lines) astrophys- 
ical S factors of 85 Rb(p,n) 85 Sr (see text). 



energy-dependence of the theoretical S factor is changed 
in such a way as to show perfect agreement with the new 
data, as seen in Fig. [3J This independently confirms the 
conclusions of previous work [5j . 

Regarding the Coulomb suppression effect, a compar- 
ison of 1.03 < / pn < 1.08 and 2.6 < / np < 3.9 shows 
that the transitions to excited states of 85 Sr are more 
important than those to states in 85 Rb in the relevant 
plasma temperature range of 2 < T < 4 GK. The almost 
negligible stellar enhancement / pn is due to the suppres- 
sion of the proton transmission coefficients to and from 
the excited states of 85 Rb for small relative proton ener- 
gies because of the Coulomb barrier. There are only few 
transitions able to contribute due to the low Q value. As 
shown by the small / pn , the transition from the ground 
state of 85 Rb dominates the proton channel. Obviously, a 
Coulomb suppression is not present in the neutron chan- 
nel. On the contrary, for this reaction / np is even more 
enhanced due to the spin structure of the available nu- 
clear levels and especially the large spin of 85 Sr ff . Be- 
cause of its large spin, it is connected to the (dominat- 
ing) low spin states in 85 Rb through higher partial waves 
than the excited states, such as the isomeric state, which 
have lower spins. Thus, the transitions from the ground 
state are suppressed by the centrifugal barrier relative to 
transitions from excited states and the latter will quickly 
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TABLE III: Astrophysical reaction rates Na{o~v) of the re- 
actions 85 Rb(p,n) 85 Sr and 85 Sr(n,p) 85 Rb computed from ex- 
perimental data. The values in italics are at temperatures 
where the experimental data mostly contribute to the rate. 
The other values are computed by supplementing theoretical 
cross sections using the modified optical potential. 
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(8.oo±U.oo)XlU 


(1. / 4±U.l 1)XIV 


0.40 
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(2.5o±U.zb J X lu 


0.60 


(1.74±0.17) x 10 


(3.49±0.35) X 10 4 


80 


(■\ 77_i_n isWin -7 


(a. sn-l-n v i n 4 


1.00 


(5.07±0.51)xl0~ 5 


(6.57±0.66)xl0 4 


1.50 


(1.28±0.13)xl0 _1 


(1.35±0.14)xl0 5 


2.00 




(2.56±0.26)xl0 5 


2.50 


(1.21±0.12)x 10 2 


(4.57±0.46)xl0 5 


3.00 


(8.22±0.82)xl0 2 


(7.81±0.78)xl0 5 


3.50 


(3.56±0.36)xl0 3 


(1.28±0.13)xl0 6 


4.00 


(1.15±0.12)xl0 4 


(2.04±0.20)xl0 6 


4.50 


(3.03±0.30)xl0 4 


(3.17±0.32)xl0 6 


5.00 


(6.89±0.69)xl0 4 


(4.76±0.48)xl0 6 


6.00 


(2.60±0.26)xl0 5 


(9.52±0.95)xl0 6 


7.00 


(7.14±0.71)xl0 5 


(1.54±0.15)xl0 7 


8.00 


(1.50±0.15)xl0 6 


(2.01±0.20)xl0 7 


9.00 


(2.50±0.25)xl0 6 


(2.18±0.22)xl0 7 


10.00 


(3.44±0.34)xl0 6 


(2.05±0.21)xl0 7 



become important, even at low temperature. As a conse- 
quence of the enhancement of / np and the suppression of 
f pn , it is more advantageous to measure the (p,n) direc- 
tion. Important transitions to states in 85 Sr are included 
in our data and the small impact of transitions from ex- 
cited states in 85 Rb is within the experimental error. 

Table [ill] gives the stellar rates for 85 Rb(p,n) 85 Sr as 
well as for 85 Sr(n,p) 85 Rb. Our data cover an energy range 
sufficient to compute the rates from 2 GK up to 4 GK. 
Because of the excellent agreement of theory with ex- 
periment, we supplement the data with the theoretical 
values to compute the rates at lower and higher temper- 
atures, applying the same errors as for the data. It is 



to be noted that fits of the rates should be obtained by 
first fitting the (n,p) rate and then deriving the (p,n) 
rate fit by modifying the fit coefficients according to de- 
tailed balance (see [3[ for details). For convenience, we 
provide the fit coefficients (including a 10% error) for 
the (n,p) rate in the widely used REACLIB format 
a Q = 16.7791+!^ ^, ax = -7.0325 x 10~ 4 , a 2 = 0.9683, 
a 3 = -9.4828, a 4 = 3.1807, a 5 = -0.3688, a 6 = 2.3328. 
The coefficients for the (p,n) rate are the same, except 



,pn _ 



17.2912 



+ ln 1.1 
+ In 0.9 



and 



-21.4342. To obtain 



the (p,n) rate the value obtained with the seven param- 
eter expression has to be multiplied by the ratio of the 
temperature-dependent partition functions. Details and 
the required partition functions can be found in Q|. 

Summary. We measured the astrophysically important 
reaction 85 Rb(p,n) 85 Sr close above the threshold in the 
energy range relevant for the 7 process. Our measure- 
ment confirms a previously derived modification of the 
global proton optical potential used in theoretical pre- 
dictions. Even more importantly, it was shown that in 
this case it is possible to derive astrophysical reaction 
rates for the (n,p) as well as the (p,n) direction from our 
(p,n) data despite of the negative reaction Q value due 
to Coulomb suppression of the stellar enhancement. A 
similar argument applies for all reactions with Q < but 
low |Q| and charged projectiles. Allowing only nucleons, 
a particles, and 7-s as projectiles or ejectiles, this effect 
still appears in more than 1200 reactions, including a 
captures relevant in the p process [H, [l(| and proton car. 
tures close to the proton dripline in the rp process |Tj 
and the vp process 
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